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Abstract 

Vibration acceleration levels on large space platforms exceed the requirements of many 
space experiments. The Microgravity Vibration Isolation Mount (MIM) was built by the Canadian 
Space Agency to attenuate these disturbances to acceptable levels, and has been operational on the 
Russian Space Station Mir since May 1996. It has demonstrated good isolation performance and 
has supported several materials science experiments. The MIM uses Lorentz (voice-coil) magnetic 
actuators to levitate and isolate payloads at the individual experiment/sub-experiment (versus rack) 
level. Payload acceleration, relative position, and relative orientation (Euler-parameter) 
measurements are fed to a state- space controller. The controller, in turn, determines the actuator 
currents needed for effective experiment isolation. This paper presents the development of an 
algebraic, state-space model of the MIM, in a form suitable for optimal controller design. The 
equations are first derived using Newton’s Second Law directly; then a second derivation (i.e., 
validation) of the same equations is provided, using Kane’s approach. 
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Introduction 


Acceleration measurements on the U.S. Space Shuttle and the Russian Mir Space Station 
show acceleration environments that are noisier than expected [1]. The acceleration environment 
on the International Space Station (ISS) likewise will not be as clean as originally anticipated; the 
ISS is unlikely to meet its microgravity requirements without the use of isolation systems [1 ], [2], 
While the quasi-static acceleration levels due to such factors as atmospheric drag, gravity gradient, 
and spacecraft rotations are on the order of several micro-g, the vibration levels above 0.01 Hz are 
likely to exceed 300 micro-g rms, with peaks typically reaching milli-g levels [3]. These 
acceleration levels are sufficient to cause significant disturbances to many experiments that have 
fluid or vapor phases, including a large class of materials science experiments [4], 

The Microgravity Vibration Isolation Mount (MIM) is designed to isolate experiments from 
the high frequency (>0.01 Hz) vibrations on the Space Shuttle, Mir, and ISS, while passing the 
quasi-static (<0.01 Hz) accelerations to the experiment [5]. It can provide up to 40 dB of 
acceleration attenuation to experiments of practically unlimited mass [6]. The acceleration- 
attenuation capability of the MIM is limited primarily by two factors; (1) the character of the 
umbilical required between the MIM base (stator) and the MIM experiment platform (flotor), and 
(2) the allowed stator-to-flotor rattlespace. A primary goal in MIM design was to isolate at the 
individual experiment, rather than entire rack, level; ideally the MIM isolates only the sensitive 
elements of an experiment. This typically results in a stator-to-flotor umbilical that can be greatly 
reduced in size and in the services it must provide, hi the current implementation, the umbilical 
provides experiments with power, and data- acquisition and control services. Even with the 
approximately 70- wire umbilical the MIM has demonstrated good isolation performance [5]. 

The first MIM unit was launched in the Priroda laboratory module which docked with Mir 
in April 1996. The system has been operational on Mir since May 1996 and has supported several 
materials science experiments. An upgraded system (MIM-2) was flown on the U.S. Space Shuttle 
on mission STS-85 in August 1997 [5], 
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Ill order to design controllers for the MIM it was necessary to develop an appropriate 
dynamic model of the system. The present paper presents an algebraic, state-space model of the 
MIM, in a form appropriate for optimal controller design. 

Problem Statement 

The dynamic modeling and microgravity vibration isolation of a tethered, one-dimensional 
experiment platform was studied extensively by Hampton [7], It was found that optimal control 
techniques could be effectively employed using a state-space system model, with relative-position, 
relative-velocity, and acceleration states. The experiment platform was assumed to be subject to 
Lorentz (voice-coil) electromagnetic actuation, and to indirect (umbilical-induced) and direct 
translational disturbances. 

The task of the research presented below was to develop a corresponding state-space model 
of the MIM. Translational and rotational relative-position, relative-velocity, and acceleration states 
were to be included, with the rotational states employing Euler parameters and their derivatives. 
The MIM dynamic model must incorporate indirect and direct translational and rotational 
disturbances. 


System Model 

A schematic of the MIM is depicted in Figure 1 . The stator, defined in reference 
frame 0 , is rigidly mounted to the orbiter. The flotor, frame 0 , is magnetically levitated above 

the stator by eight Lorentz actuators (two shown), each consisting of a flat racetrack- shaped 
electrical coil positioned between a set of Nd-Fe-Bo supermagnets. The coils and the 
supermagnets are fixed to the stator and flotor, respectively. Control currents passing through the 
coils interact with their respective supermagnet flux fields to produce control forces used for flotor 
isolation and disturbance attenuation [8]. 

^ zv z\ 

The flotor has mass center F and a dextral coordinate system with unit vectors / , / , and 
/ , and origin F 0 . The stator (actually, stator-plus-orbiter) has mass center S and a dextral 
coordinate system with unit vectors $, s , and $ , and origin S 0 . The inertial reference 
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frame © is similarly defined by fi , h , and h , and origin A/ 0 . The umbilical is attached to the 

stator at S u , and to the flotor at F u . When the flotor is centered in its rattlespace (the “home” 
position), F and F u are located at stator-fixed points F * , and F ih , respectively. 



Figure 1 . Schematic of the MEM 


State Equations of Motion 

Translational Equations of Motion 

Let E be some flotor-fixed point of interest for which the acceleration is to be determined. 


If E has inertial position r , then its inertial velocity and acceleration are f = — r and 

r -N 0 E J -N E (If \-N E/ 


d\ N d 


~ N E d t d t 


r N , respectively. (The presuperscript indicates the reference frame of the 


differentiations. The subscripts indicate the vector origin and terminus.) The angular velocity and 
angular acceleration of the flotor with respect to the inertial frame are represented by N of and 


V , respectively, where N a F = ■ 
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Let F be the resultant of all external forces acting on the flotor; M_ F 1 (or simply M), the 
moment resultant of these forces about F\ m, the flotor mass; and 1 (or / ), the central inertia 


dyadic of the flotor for / / , and f . Then Newton’s Second Law for the flotor can be 


expressed in the following two forms: 


F = rrif , 


and 


M_= F t h F + %/x(/ % F ), 


(Eq. 1) 
(Eq. 2) 


From Equation (2), 


t h F = r l \M-% F x(i% F ). 


(Eq. 3) 


It will be useful to find an expression for r_ N E in terms of the acceleration r_ N s of the umbilical 
attachment point S u , and in terms of the extension of the umbilical from its relaxed position. 


Begin with the following: r, F =r_ NE + r_ EF - r N So - r, o . 


(Eq. 4) 


Differentiation of Equation (4) yields 


ts F E + ® X Lef i N S 0 — X -S 0 S 


(Eq. 5) 


A second differentiation gives 

: L No e+ V xr ff + ¥x( N co F xl ef )-l ns -Yxlss- Y 


r 

-S F 


S.A N S 

— "i-is / • 


x| "©“xr, c . ) . (Eq. 6) 


Substitution for N a F from Equation (3) into Equation (6) yields 




M- N ( 0 F x 


- N co s x N co s xr ss + N co F x N co F xr FF 

— ono — thr 


(/• A m F | Jxr £F -l NS - 

( 


Y xr ss 


(Eq. 7) 


hi these equations 


N ^ I NS . F 

0) = 0) + 0) . 


(Eq. 8) 
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Under the assumptions that !V ffl s and N a s are negligibly small and, therefore, that 


r ~ r 
—N S —N S u ’ 


(Eq. 9) 


Equation (7) reduces to 


y = y + 1 / 1 • 

-S„F„ -N„E ^ > 1 


M_- s 0 ) F x (/■ W'jjj 


* r -EF £-NqS — * ( — * - EF 


(Eq. 10) 


Linearization about s O) F = 0 yields the following result: 

Ls F — iijv E {L 'M^\^LeF — Ln s ■ (Eq. 11) 

Appropriate expressions for F and M will now be detenuined, for substitution into Equations (1) 
and (11), respectively. Those equations will be used in turn to obtain a more useful expression for 
y s F . (See Equations (43-48).] 


I 

The force resultant F is the vector sum of the eight actuator (coil) forces F (/ = 1, ...,8), 

with resultant F c ; of the umbilical force F„,, caused by umbilical extensions from the relaxed 
position; of the direct disturbance forces, with resultant F d ; and of the gravitational force F g - 
Gravity may be neglected for a space vehicle in free-fall orbit. The moment resultant M is the 
vector sum of the moments due to the coil forces, with resultant M ; of the moment M , due to 

the umbilical force F ■ of the moment M due to the umbilical rotations from the relaxed 

— ut 7 ur 

orientation; and of the moment M_ d due to the direct disturbance forces. There is no moment due 

to gravity, since M is about the flotor center of mass F . hi equation form, assuming the i th coil 
force to be applied at flotor-fixed point B,, 

1=1 

8 

and K = b x E.' c +L f F xF, <t +K lir +K d - (Eq. 13) 

i=i 
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More explicit expressions for F_' c and F_ ut will now be developed. If the actuator has coil 
current /, / ; , length L. , and magnetic flux density B j B j , then the associated actuator force 
becomes 


F\ =-/ I. B / , x />\ . 


(Eq. 14) 


Assume a translational stiffness K' t for an umbilical elongation in the s. direction, and a 

corresponding translational damping C \ . Let F_ b represent the umbilical bias force, exerted by the 
umbilical on the flotor in the home position. Then the total force of the umbilical on the flotor 
becomes 


F ut = -\ f ~Lsf )A,k + X C < T ( Ls F ~ Ls F ) ' ^ 

1=1 L /= 1 c 't 


I, > + F h 


(Eq. 15) 


Define the following, for i = 1, 2, 3 : 


x ai = (Ls f — Ls f )•!; (Eq. 16) 

and x bi = x aj . (Eq. 17) 

If V =0, Equation (15) becomes 






+ F„ 


(Eq. 18) 


L '=i J 

The relative positions x ai and the relative velocities x bi will be six of the nine translational states 
used in the state-space formulation of the system equations of motion. 

As with F_' c and F_ M above, M ^ can also be expressed in more explicit form, in analogous 

fashion. Assume a rotational stiffness K' r and a rotational damping C' , for umbilical twist about 
the s direction. Let §'' ls n represent the rotation of the flotor, relative to the stator, from the 
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relative position in which the / and s coordinate systems are aligned. F /^is the rotation axis, 

and 4> is the angle of twist about that axis. Let M ,, represent the umbilical bias moment, exerted by 
the umbilical on the flotor in the home position. Then the moment M_ ur can be expressed by the 
following: 



(Eq. 19) 
(Eq. 20) 


Equation (20) can be expressed in alternate form using Euler parameters. Let F/S n be 
described in © by 


F '% = + + ^ 3—3 ■ (Eq. 2 1 ) 

Define the following Euler parameters [9], 

F/S p 0 =cos|, (Eq. 22) 

[S)F %= ei sm^, (Eq. 23) 

(S)F/5 P 2 = e 2 sin|, (Eq. 24) 

(5)F/5 p 3 =e 3 sin|, (Eq. 25) 

and F/5 |3=sin^ F/5 n • (Eq. 26) 

— Z —(j) 

For small values of (|), the Euler parameters can be simplified: 

F/s Po=l> (Eq. 27) 
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(W P, =e^/2, 

™%=eA/2, 


and 


Note that, for small angles. 


Fls V=l F/s h ■ 

— £ ~<P 


V%=2 F % 

This equation can be used to simplify the stiffness terms of Equation (20). 

As for the damping tenn, Equation (31) can be differentiated to yield 


F / S - ^ ___ ^ 


R = —cos A F/S n + sin A f/sx 
” ^ ~ _0 O 'A ’ 


. 


or, for small angles, 


Equation (20) now becomes 




. F / Sa 


2 ‘ p=0' V/ fd' n, 
— — 0 1 — 0 


M .... = -2 


1 = 1 i=l ' ' 


+ M , 


Define the following, for / = 1, 2, 3 : 


- =' 'I* i 


and 




The assumption that N or is negligible yields, finally. 


^=-2 \t[Kx* + CXM+M t 


i=i 


(Eq. 28) 
(Eq. 29) 
(Eq. 30) 

(Eq. 31) 


(Eq. 32) 


(Eq. 33) 
(Eq. 34) 

(Eq. 35) 

(Eq. 36) 
(Eq. 37) 

(Eq. 38) 
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Note that Equation (11) describes r‘ ,, in terms of the acceleration of an arbitrary flotor- 

u*u 


fixed pomt E. For E located at flotor mass center F * Equation (11) can be used straightforwardly 
with Equation (1) to yield 

Ls F = J t E- + 1" 1 ■ M) x r - F ' F ~ L No s • (Eq. 39) 

Define now three unknown-acceleration tenns, to be used with Equation (39). The first term 
represents the indirect translational acceleration disturbance input to the flotor, applied at the stator 
end of the umbilical: 


«/„ = 


-N 0 S, t 


(Eq. 40) 


The second term represents the direct translational acceleration disturbance to the flotor, due to 
unknown disturbance force F_ d : 


a d =~F d 

m 


(Eq. 41) 


And the third represents the direct angular acceleration disturbance input to the flotor, due to F_ d 


Qd~ L ' — d ’ 


(Eq. 42) 


Substitution from Equations (12), (13), (14), (18), (38), (40), (41), and (42) into (39) yields the 
following result: 

is uFu = ~^c+ —ui ) + L~ 1 - (M c + K ul + K ur )xr FF + a d xr 


-F F —in T — d ' 


8 8 


where 


£. = XE = E(-wu<x£), 


/=! i=l 




t( K X,+cx,)s, 


+ E_ h , 


«, = 'Em!, = Xz>»,x(-/,w, xb,), 

i = i ;=i 


(Eq. 43) 


(Eq. 44) 


(Eq. 45) 


(Eq. 46) 
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(Eq. 47) 


and 


M-m ~ X 


S(^>.,+c>„)s,. 

1=1 


M ur =-2 £[*'**+ C^]*, \+M b . 


(Eq. 48) 


L i=i J 

Substitution from Equation (43) into Equation (11) produces the following equation for the 
acceleration of arbitrary flotor point E: 


r =—(f +F ,) + r 1 -(m +M , + M )xr„ + 

Assuming N (O s to be negligible, one also has the following: 

s d 

Lsf ~ dt 


QLd -^Lfe + (Ld 


(Eq. 49) 


and 


s d 2 

E s f ~ df 2 


(-S F ) ’ 

(Eq. 50) 

)‘ 

(Eq. 51) 


(Note that assuming N (o_ s to be negligible does not imply that © and © are identical; it means 

rather that © can be treated as if it is in pure translation relative to © for the frequencies of 

interest.) Equations (43), (49), (50), and (51) provide the basis for a state-space form of the 
translational equations of motion, using x ai , x bj , and low-pass-filtered approximations to the 

components of r_ NaE [see Equations (94) and (99)], as states. 


Rotational Equations of Motion 

Let ^(3= sin — F/S n as before [Eq. (26)]. Differentiating the left side twice produces 

— 2 dj 


d 

dt 


( ,/A p)= ( ( l H) + A m v x ,/s j3 


(Eq. 52) 


N , 2 


S . 2 


and ^r( F/5 p) = ^r( F/5 p)+V x-^-(^p)+% 5 x F/s p + V x F/S p. (Eq. 53) 


dt 


dt 
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Assuming as before that N a s ~ 0 , Equations (52) and (53) become, respectively, 
and F ' S P=-^r( F 'f)- 


Returning to Equation (26), two differentiations of the right side yield 



U.Q+ 4>cos^- n^+ smj 


(Eq. 54) 
(Eq. 55) 

(Eq. 56) 


Linearizing about <f> = 0 and <j> = 0 , Equation (56) becomes 

2 = $ F,S n , + 2 <j> F %+ b F % . (Eq. 57) 

Equations (34), (54), (55), and (57) provide the basis for a state-space form of the rotational 
equations of motion, using as states the s. components of 7 s |3 and of [) (i.e., x di and x ei , 
respectively, for i = 1, 2, 3). 


Equations of Motion in State-Space Form 
From Equation (16), 

Ls F ~Ls F = + X a 2§-2 + -*a3^3 • (Eq. 58) 

Differentiation, along with the use of Equations (17) and (50), leads to the following: 

Ls F — Ls F = ^ b\ — 1 + ^b 2^2 + ^b3^-3 • (Eq. 59) 

A second differentiation yields 

Ls F ~ Ls F 1 + -* 62^2 + A,3^3 • (Eq. 60) 

Introduce the use of a pre superscript in parentheses to indicate the coordinate system used for 
componentiation. (This notation allows vectors to be expressed unambiguously in terms of their 
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measure numbers.) Then Equations (58) and (59) take the respective fonns, 


x. 


a 1 


(S) 


(S) 


x 


a 2 


= x. 


(Eq. 61) 


x 


a 3 j 


and 


(S). (S). (S). 

j* m __ v m ^ y y 

Ls F Ls F — Ls F — ±a — ±b’ 


(Eq. 62) 

where x a and x b are defined as indicated. x d and x e have corresponding definitions. [Cf. 
Equations (36) and (37).] 

Equations (43) and (60) can be used together to develop a state-space equation for x_ b . 
First, express Equation (60) in measure-number fonn: 




—S U F„ ~ ' 


X 


bl 


-V, 


bl 


X, 


bl 




(Eq. 63) 


Next, define rotation matrix S F Q by 


h 


f > 

/V 

L 

/V 

< v 

—2 

■ = 5/F e ' 

/V 

L 

UsJ 


/V 

[L] 


(Eq. 64) 


where the prefix indicates the rotation of frame © relative to frame © . Finally, observe that, for 


arbitrary vectors 

Li = •e/ 1 + yj 2 + Zl / 3 
and r 2 = x 2 f^ + y 2 f_ 2 + z 2 / 3 , 

the cross product can be expressed in determinant form by 




L 

L 


U XLi = 

X 

yi 



*2 

y 2 
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(Eq. 65) 
(Eq. 66) 


(Eq. 67) 



or in matrix form (i.e., using measure numbers) [9], by 


(F) 


fcl Xf 2 ): 


‘ 0 

-A 

Vi 


x 2 

Z 1 

0 

--E 

< 

V 2 

_-A 

-*) 

0 _ 


. Z 2. 


(Eq. 68) 


Represent the above skew-syiumetric matrix by [F> r * . Using this notation, Equation (43) can be 
expressed as follows: 

1_ s/Fn H mWaa ]- s l F n _<«/» 


* b = ~l ±U+ <2 "A£+^+^|-^ Iff/ M d ~ Q in + K a d , (Eq. 69) 


where 7 is the inertia matrix corresponding to / . 


Linearizing Equation (3) about N co F = 0 yields 

N a F = I~ l -M. 


But 


N F 

a =■ 


Njl 

dr 2 


(o ' S Z7 0 ) =0 7 '/> 0 + 2o ? '% + 0 ; '/?„ 


From Equations (57), (70), and (71), 


or, equivalently, 


2 F % = r 1 K, 


V 

dr 


( F/ 1P) 


r l ■ m , 


Application of Equation (55) leads directly to 

dr v 2= — 

hi measure-number form, 

L=\ s/F Qr l(F) K, 

or, equivalently, i, = | 5/F 0 /“'[ (F) M ut + [F) M iir + (F) M c ] + ^ S "'Q {F h d . 


(Eq. 70) 
(Eq. 71) 

(Eq. 72) 
(Eq. 73) 

(Eq. 74) 

(Eq. 75) 
(Eq. 76) 
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Six state equations of the system are given by Equations (17) and (37), iterating on /; six 
more, by Equations (69) and (76). The latter six are written in terms of the various forces and 
moments acting on the system, which loads have been defined in vector form by Equations (44) 
through (48). These loads can be rewritten in measure-number form and substituted into Equations 
(69) and (76), as follows. Beginning with Equation (44), the i ,h control force can be expressed as 


(S)pi 


(S)~x (F)~ 

-Z, /, s/f Q B t B, 


h =FcUi, 


(Eq. 77) 


The resultant control force becomes 


(5) 


] F c =Y J (S) F i c = Fu 

C X—md C C ’ 


(Eq. 78) 


i=i 


where F ‘ , F c , u, , and u are defined as indicated. 


Next, using Equation (61) with (45), the translational force the umbilical exerts on the 
flotor can be expressed by 


(S) ■ 


Lu, = -K,x a - C,x b + ' S) F b = F uta x a + F utb x b + '*'F b , 


(s > , 


where 


K, 


K) 0 0 

o k ; o 

0 0 K , 3 


(Eq. 79) 


(Eq. 80) 


C, 


c! o 


o 

o 


c 


0 

0 


0 c. 


(Eq. 81) 


and F uta and F mh are defined as indicated. 


The i th control force F_‘ c exerts on the flotor a moment Af , defined by Equation (46). 


Using again the notation introduced with Equation (68), this moment can be expressed by 

I /, = M; Ui , (Eq. 82) 


{f) M[ = 




i C l 
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and the resultant moment by 


] M = t {F) M i = M u, 

c c c — ’ 


where M ' and M are defined as indicated. 


The umbilical force F_ ut exerts on the flotor a moment M , given by Equation (47). 
Substituting from Equation (79), this moment can be expressed by 


{F) M = (F) r x . S I F O 1 \ F x +F x 1' 

JlLut -F F„ \L V r ulaXLa ^ ‘mb —b J ‘ 


(Eq. 84) 


or, alternatively, 


Mu, = M uta x a + M utb x b , 


(Eq. 85) 


for M uta and M utb appropriately defined. 

Finally, Equation (48) expresses the moment M :r that the umbilical applies to the flotor due 
to umbilical rotational stiffness. The following equations express M , r . in measure-number fonn: 


] M W . =-2 s/F Q T K r x d -2 s/F Q T C,.x e+ [F) M b = M urd x d + m M b 


(Eq. 86) 


where 


'K\ 0 O' 

OK; 0 , 
0 0 Kl 

'C\ 0 o' 

0 C r 2 0 , 

0 0 cl 


(Eq. 87) 


(Eq. 88) 


and M urd and M ure are appropriately defined. 

Substituting from Equations (77) through (88), Equations (69) and (76) become, respectively, 


.p - S / F n (F) r x I~ l M Iv + — F - s/f O {F) r x I~ l M Iv 

r uta Xd LFF 1 m ula tla ^ r mb Xd Lf F 1 lV1 utb 


+(- s,F Q ' F y FF r‘M lM )x J+ (- s/F Q m i rr r'Mj)x,+{-F r slr Q 'Y ff 1 


M, \u 


(Eq. 89) 


+-®&- s ' f e ,r y„ ' F 'q J - SIF Q m y. r r' m M„ 
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and 


le -T; SF Q 1 i [ M u,a^a + M utb*b + M urd*d + M ure*e] 

+ L SIF Qr l M c u+^ s/F Q (F h d +^ s/F QI~ l {F) M b . 

For completeness, Equation (37) can be rewritten as x d = x e . 

To include (s ’l No i; as states, define x c by 

a h S) L NoE =L + (O h x c , 

for some high value of circular frequency Ct) h . Taking the Laplace Transform, 



so that x c ~ (S ] r_ NoE for co« 0 ) h . 


(Eq. 90) 
(Eq. 91) 

(Eq. 92) 

(Eq. 93) 
(Eq. 94) 


Now using Equations (78), (79), 83), (85), (86), and (92) with (49), 


x - co f— F - S/F 0 (F V X . lx + a> {— F - s,f O <F) r x , / lx -to x 

±c — Lu h r um v: LFE 1 m uta ±a +UJ h r utb ^ '~F'E l LV1 utb ti b W h±c 

\m ) V m r 


-co h ( s,F Q (F) i FE r l M w . d )x d -to h ( s/F Q {F y FE r l M ure )x e +to ,[^ f- sif q {F V FE r l M c y 95 ) 
+to h ( % +©/-l <S) £ i -to„( SIF Q iF Y FE ) (F) cc d -co h ( SIF Q (P) r x PE r l ) (F) M b . 


A state-space representation of the system is given by Equations (62), (89), (90), (91), and 
(95), for state vector 



(Eq. 96) 
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For the small rotation angles associated with the MIM, s r Q is approximately equal to the 3x3 
identity matrix, in which case the state equations have constant coefficients. Specifically, 


X.a=* b ’ 


- F- {F) rl j 7 r l M,„„ V + { - F„„- (F) rl , r l M.„ 


1 uta —F F uta t! a 1 

m j l m 


-FF 1,1 lUh lib 


{ m £ r r'Mj l x d -( ]x, +(-F c - m r* rF r' 


M r , \u 


(S) _L (S k, _l_ ^ (S) Z7 (F > X (F> SV (F) X 7-1 A 

- Q.i„+ Q d + — E b ~ r FF a d - r FF I M b , 
m 


i =oJ-F„- A* £ /->m_ h k -<o„x c 

\m r v m r 


-®4 hj A'M-k +®. - F - 

\ / \ ’ \m 

*», •% +4- L >’s -®; F, £k(’ F1 a +/-' '"M.) ■ 

m \ > 


Xd = x e , 


^e=-^ l [ M , lla -±„ + M ulb Xb + M uvd*d + M are* e\ 

+\r l M c u+\^h 1+ /- mf ’m s ). 


(Eq. 97) 
(Eq. 98) 

(Eq. 99) 

(Eq. 100) 
(Eq. 101) 


State-Space Equations With MIM-2 States 

The above state equations must be modified to account for the states actually used in 
MIM-2. Designate by the post-superscript R the states defined above; and by the post-superscript 
C the states used with MIM-2. The fonner set of states are as follows, for i = 1, 2, 3 : 

x R ai =L Fuh F u -ii, (Eq. 102) 

xi = x« , (Eq. 103) 
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1 


(Eq. 104) 


and 

and the latter, 





(Eq. 105) 
(Eq. 106) 
(Eq. 107) 
(Eq. 108) 

(Eq. 109) 

(Eq. 110) 
(Eq. Ill) 

(Eq. 112) 


Differentiating twice, under the prior assumption that N (O s is negligible, leads to the following two 


equations: 

dL-S F d Hf'f ' F 

— + O) X [_f'f 

dt dt 

(Eq. 113) 

otin 

■5 j 2 s H 2 r » » 

a Ls u F u U -F h F s F S F j S F \ 

(Eq. 114) 

allU 

2 — 2 + X ’_i f + Us X 1 ^ i_/. f J • 

tt i tc r 

Linearizing about 

n g/ = 0 as before, the following interrelationships are found to hold for the two 

sets of states: 

x:=xU( S/F Q (F) Ly F+ ( %y), 

(Eq. 115) 


d=£~2 S/F Q (F) rJ^ lF Q r x« , 

(Eq. 116) 


£ = £- SIF Q m i> F r UF) M. 

(Eq. 117) 
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Designate by ! ' ) r b the final two terms of Equation (115). Substituting now from Equations (115), 
(116), and (117) into Equations (62), (89), (90), (91), and (95) yields the following state equations 
for MIM-2: 


. C C 

I a =l b > 


(Eq. 118) 


• c 

X, = 


If ] 

£+ 

[If ) 

£+ 

\m V 


J 

V 




x . 


+ 


-I 

\m ) 


M- ( V+ (S, fl,+ 1 


'1 F ' 

, uta , 

\m J 


(S) 


r b + 


\m) 


F 

L-b > 


(Eq. 119) 


x c =(o\-F - s ' F 0 (F) r x I~ l M V+flil— F - 5,f O (F) r x F l M 

~ c lU \m ruta ^ ~ FEjL Lluta j m Futb ^ ~ F r1 1 


.c 

mb \*-b 


-C0 h x c c+ C0 h {- s l F Q {F Yr E r l Mjx 


-2 


+«* — Fj F Q (F) l ff „ sif q t +2 s > f q < f Y fe YmJ f q (F Y FFu s,F Q T - s,F Q (F Y F - E r l M, 


m 

«>4 Y~ !IFq 


c 

we ttle 


(Eq. 120) 


+m {i F -- VQ ' F V‘ r ' M Y'-‘ +m F"Q ,F %A m & + r‘ <F> Mb) ■ 


■ c c 

Id = le ’ 


(Eq. 121) 


. c 

x. = 


+ 


+ 


Y s/ ~‘ m " 


V/ 


j 


— s/f qi~ 1 m 
2 


utb 


lb + 


- SIF Qr'M^ 


Id 


\ 


SIF Qr y MjF Q «Y f . f s,f Q t 


J 


i e + | \ s ' F Qr l M c 


6 1 


s "QF 1 m u 


(S) 


Lt + 


(1 


S/F Q [ F k d + r UF) M b ). 


(Eq. 122) 


Again, for the small rotation angles associated with the MIM, s F Q is approximately equal to the 


3x3 identity matrix; Equations (118) through (122) reduce to the following: 


. C C 

la =I b > 


(Eq. 123) 
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a = 


'i p v , 

— r V 4 - 

uta la 

\ m J 


V ( 2 


~ F utb k + 

m 


- (S) a + (s> a ■ 

'±itT r H rf 


{If ) 


+ 

rr 

[m J 





— F (F) r x , 

J mb LF F u 

m 


|(S) Z 


V (1 ' 


V + 


-F f U 

\ m J 


(Eq. 124) 


■ c 

X ' = CO, 


(F)X r-1 


A 


Lr ,I M u ,„ k, +C0 


m 


uta LF E uta 


m 


T - (F) x r-l,, 1C „ C 

" F utb Lf'E^ M utb k —CO h X r 


+co 


U F Y, i~ ] m V c 

i\ LF E* Lvl urd}±d 


+C0, 


— F ut b (F Yf>f + 2 (F Vk - <F) rYr l M.. 


m 


+co. 


+ 0 ) 


-f c - ( YfY~ 1 k\w%+^ 

m J- 

l-Y E \ (F) a d + r^M b ), 


LF'E 1 lr± utb LF'F ■ I-F'E* ,u .k tie 


— F - <F) r x . 7 _1 M PV 

1 uta Lf E 1 1V1 uta Lb 
m 


pi 



V 


(Eq. 125) 


.c c 

= I* , 


(Eq. 126) 


• c 

— e : 


+ 


— 7 _1 M 
2 


J— fl 


n£ + 




& + 

V 2 


2 y 


-r l M...-r l M... u (F V X , 


\ 

( 1 3 

-r l M c 

x c e + 

) 

2 

V ^ / 




(S) 


D>+ 


II< F '« i+ /- ,<F v s ). 


(Eq. 127) 


State-Space Equations Using Kane’s Dynamics [10] 

The above state equations for MIM-2 can be derived alternatively, by the approach 
commonly called “Kane’s Dynamics.” First define generalized coordinates q, and generalized 
speeds as follows. [Note the use of post-superscripts now, instead of the previous post- 
subscripts, with position (and, later, velocity and acceleration) vectors. The post-subscript 
position, with vectors, is used for another purpose in Kane’s notation, as will be seen.] 


<h =L 
Qi = L 

<?3 =L 


Kf’ 

F' h F' 

F' h F’ 


■ s 


1 ’ 


-2 ’ 


(Eq. 128) 
(Eq. 129) 
(Eq. 130) 
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(Eq. 131) 


N F * 

u 4 = ® s 1 =q 4 , 

u 5 = N u F -l 2 =q 5 , (Eq. 132) 

u 6 =V-s 3 =q 6 , (Eq. 133) 

m 7 =<?? =- c ° h qi+° ) hk NoE -ii)’ ( E q- 134 ) 

= q, = -Oi h q^ + 0 ) h {y N ° E • s 2 ), (Eq. 135) 

u 9 =q 9 = ~0) h q 9 + (O h (y N,)l ' : • s 3 ), (Eq. 136) 

u w= N ®. s ■ li = < 710 ’ (Eq- 137 ) 

u n = N (O s -s 2 =q n , (Eq. 138) 

u i 2 = N (O s -i 3 =q l2 , (Eq. 139) 

Mi = (h > (Eq- 140) 

u 2 =q 2 , (Eq. 141) 

and u 3 =q 3 . (Eq. 142) 


Next determine useful velocities and angular velocities in terms of these generalized 
coordinates and generalized speeds. The angular velocities can be expressed as follows: 

N 0) s = m ^ Ij + u u s 2 + u i2$i ’ (Eq. 143 ) 

N qf = m 4 1 1 + u 5 s _ 2 + u 6 s _ 3 , (Eq. 144 ) 

and s a/ = N af - 'V = (n 4 - m l0 )| 1 + (n 5 - m 11 )| 2 + ( m 6 - n 12 )| 3 . (Eq. 145 ) 

The velocity of the flotor center-of-mass is 


v F = 


V 

(It 


,.N„S 


. S F . F F 

+ r + r 


S . AT S w S F 

- v + a) x r 


r/r 


/ F F \ . W 

lr 1+ m 


S w F F 

x r 


(Eq. 146) 
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Iii terms of the generalized coordinates and generalized speeds, 


/ = v s + W x r s F + y, ( u x + r/ 3 , - q 2 u vl ) 


^2V^2 #1^12 #3^10 


) + | 3 (m 3 


+ <7 2 w io <7i% 


(Eq. 147) 


Accordingly, the respective velocities of the isolation point E , the eight actuator force-application 
points B I , and the umbilical-attachment point F u , are as follows: 


E S N..S S F ‘ ~ l 

v =v + a) xr 


+ ^3% -<72%) 

+ 5 2 (^2 <7l% _ <?3%) <?2«10 _ <7l%) 


+ N af xr f£ , 


(Eq. 148) 


and 


v B ‘ =v‘ s “+ W xr A ' f ‘ + i (zq + r/ 3 M u -q 2 u ii) 

+ 5 2 (m 2 + <?!% — <73%)"^ ^3 (^3 + <7 2 M 10 _ <7l%) 


N F F B 

+ cq xr \ 


y = y s " + x (iq + r/ 3 w u - < 7 2 zq 2 ) 

+ 5 2 (m 2 + <7l«l 2 — <?3% )e £3 (« 3 + <?2% — <7l% ) 


, N F F F u 

+ Q) xr “ 


(Eq. 149) 


(Eq. 150) 


One can now express the linearized partial velocities (L.P.V.’s) and linearized partial 
angular velocities (L.P.A.V.’s), corresponding to the foregoing velocities and angular velocities, 


using the following notation: yf 


dv E 


is the partial velocity (P.V.) of point E with respect to the 


i ,h generalized speed u, , N 0 ) F = — - — is the partial angular velocity (P.A.V.) of reference frame F 

a Uj 

with respect to reference frame N, and ,vf and fmf are the respective linearized velocity tenns. 
Assume now that v 0 ) s ~ 0 , so that w 10 = u u = u V2 = 0 . Then yf" = 0 (/ = 10, 1 1, 12) . Assume 
additionally that the flotor mass is much smaller than that of the combined stator-plus-orbiter, in 
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which case s O) 1 does not affect v s, ‘ . Then vf“ = 0 ( / = 4,5,6) . Similarly, since y' s " is not affected 

by r F,y or r N ° E , vf- = 0 (i = 1, 2, 3, 7, 8, 9) . That is, 

vf- =0(/ = l,...,12). (Eq. 151) 

Fnrther, since y 5 " is unaffected by r_ F E , r F B ‘ , or r F ! “ , one has the following L.P.V.’s and 
L.P.A.V.’s: 


,vf =li (i = 1, 2,3), 

(Eq. 152) 

,vf =0 (/ = 4,. ..,9), 

(Eq. 153) 

vf =/ vf ' =,yf“ =,vf V; in {l,...,9}and V/ in {l,...,8}, 

(Eq. 154) 

ftof =0 (z = l, 2, 3, 7, 8, 9), 

(Eq. 155) 

/ ^4 - li ’ 

(Eq. 156) 

S..F £ 

/C() 5 — 5 2 > 

(Eq. 157) 

S..F j 

/^6 — ^.3 ’ 

(Eq. 158) 

AF . F S_F /• 1 

= ,&>, 9). 

(Eq. 159) 


All L.P.V.’s and L.P.A.V.’s associated with u w , u n , and u l2 are 0. 

hi order to determine the linearized accelerations (L.A.’s) and linearized angular 
accelerations (L.A.A.’s), one must first determine the linearized velocities (L.V.’s) and the 
linearized angular velocities (L.A.V.’s). Assuming still that V 0 ) s ~ 0 , the L.V.’s and L.A.V.’s are 
as follows: 

,v F =v_ Su +u 1 s_ l + u 2 s_ 2 +u 3 s 3 , (Eq. 160) 

E S u a /v a j 5 F F F | /~i — « 4 s- a \ 

,v =v +u l s l + u 2 s 2 +u 3 s 3 + l \(0 xr I, (Eq. 161) 
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B i S u , x , X x i S F w I 

,v =v +u 1 s l +u 2 s 2 +u 3 s 3 +i\ 6 ) xr ), 


(Eq. 162) 


F„ s„ , X „ X / S F F'F U 1 

; v =v +u x s { +u 2 s 2 + u 3 s 3 +,y (Q xr J, 


(Eq. 163) 


N, n 

, m =0, 


(Eq. 164) 

and 

N F S F A A A 

,CO = ,(0 = U 4 S l + U 5 l 2 + u 6 s 3 . 


(Eq. 165) 


Differentiating the L.V.’s and L.A.V.’s, and linearizing about S (Q F 

~ 0 , yields the 

following L.A.’s and L.A.A.’s: 




,a F = a s “ + + ii 2 s_ 2 + M3I3, 


(Eq. 166) 


E S a yv a./a. a. ^ > 

jU = a “ + + u 2 s_ 2 + it 3 $ 3 + (w^ + u s S 2 + h 6 i 3; 

|x/ £ , 

(Eq. 167) 


,a B ‘ = + w 2 s 2 + h 3 £ 3 + (m 4 s 1 + m 5 | 2 + i/ 6 | 3 ] 

I x r FBi 

(Eq. 168) 


F A A A, a/a, A a ) 

" = a " + hqlj + h 2 | 2 + h 3 l 3 + + m 5 S 2 + m 6 5 3; 

)xr FF “ 

(Eq. 169) 


Fas* 

C/> 

II 

IO 


(Eq. 170) 

and 

N F S F .A .A .A 

,a =,a = w 4 ^i + W 5 I 2 + m 6^3- 


(Eq. 171) 


Beginning with Equation (3), one can also obtain readily an alternate expression for fa' : 

x(r„-rj]+(li:" xf;| (Eq. 172 ) 

The final step, before writing the generalized active forces and generalized inertia forces of 
Kane’s equations, is to determine the contributing loads: 

The resultant of the actuator forces (cf. Eq. 44), which are considered to be applied at respective 
locations B i , is 

F, (Eq. 173) 

/= 1 i=l 
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The umbilical force (cf. Eq. 15), with the former term (in curly brackets) applied at F u , and the 
latter (umbilical bias force F _ b ) at F* , is 


1=1 


\- K ‘(r F ^-s\ 

—C‘ 

-i '-1 

"n, 

- (r FukF “ )• s j 

• 1 
s i 


clt ’ 

J 


■+£*■ 


(Eq. 174) 


The direct disturbance force (cf. Eq. 41), applied at F , is F_ d (unknown). 
The umbilical moment (cf. Eq. 20), applied about F* , is 

Kir = t {-<[(< l ' , %) i ] i ,- c ' r [ h ,/ s 2 ,)- s , 

i=i L 


l +K b . 


(Eq. 175) 
(Eq. 176) 


The direct disturbance moment (cf. Eq. 42), applied about F , is M_ d (unknown). 

The above expressions for the contributing loads can also be written in measure-number 
fonn, needed for eventual computer implementation. Equation (173) becomes 


(5) 


F = 


-L^T^QBfX 


i 1 + ... + \-l 8 {S T s s/f qb s ( % 


where / = [/ 1 ,...,/ 8 ] r . 

Equation (174) is first re-expressed as 


h=F c L, (Eq. 177a, b) 
(Eq. 178) 


Em - + / F - (it -c;[”4 ^ 

i=i dt 


r “ A +r FhF +r FF “ 1-s , 


l > + F h 


(Eq. 179) 


This now can readily be rewritten as 



Qi 


M 1 


w 4 

(5) f = p 

— ut uto 

Qi 

' + F 

utb 

u 2 

,+ <S) p ^ _ p S /I'Q F'q! - 

m 5 


( h_ 


u 3 


m 6 . 


t + F„ 


( s ) 


r b . (Eq. 180) 


The direct disturbance force is simply lS) F d . 
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Equation (175) becomes 


(F) 


M... = M u J 


Qa 


u 4 

Qs 

■+M urr / 

u 5 

q 6 . 


u 6 


(■ F)t 


where M urd = - SIF Q T K r 

and = - s,F Q T C r . 

The direct disturbance moment is M_ d = I / • a A . 


(Eq. 181) 

(Eq. 182) 
(Eq. 183) 

(Eq. 184) 


Using the expressions for the L.P.V.’s and the L.P.A.V.’s with the contributing loads, the 
generalized active forces Q i can be determined as follows: 


a 


,vf •(£,+£») 


,vf- •(£. -£») 

8 

'/vf -E~ 

a 

> = < 

, if •(£„+£») 


•(£. -£„) 

+s 


.a. 


, if •(£,+£») 


, vf ■(£,,-£.) 

i=l 



- [S ^p + [S) F + (S) F 

— — ut ' —€■ —d ? 


a 

a 

a 


k:-(K ur +M d+ 


K ■{K ur +K d + 


l®6 • ( M-ur + K d + 


r F F x 


F F 

r x 


r F F x 


(£„-£.)]+ fie" xf; 

V i=l 

(£.,-£>)]+ fix" xf; 


(Eq. 185a, b) 


(S) 


F "x(F„ 



'xF'' 


+ 


s,f q( 


(F) 


M 


{F) m 


\ 




a 


0 

and 

■< 

a 

> = < 

0 



a. 


0 


(Eq. 186a, b) 


(Eq. 187) 
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The generalized inertia forces are 


q: 


',vf 

(- m l a F ) 



ft* 

> = - 

,v£ 

(- m l a F j 

■ = -m (s y\ 

(Eq. 188a, b) 

.ft*. 



(- rn l a_ F j 




r 

ft 




_ * 

ft 

► = < 

.,W) 

_ * 

.ft. 


K-(-I-fa F -f© F x/-,V) 


= sif Q {- 1 ■ S o L ' ) (linearizing about =0), (Eq. 189a, b) 




tel 


0 

and 

< 


► = < 

0 



.ft*. 


0 


(Eq. 190) 


Kane’s Dynamical Equations for MIM can now be written in the following, recognizable 
fonns of Newton’s Second Law: 


(S) L ul + (S) L c + [S) L d -m iS) l a F =0, 


(Eq. 191) 


where [cf. Eq. (166)] 


(s la F = (S) a in +i 


(Eq. 192) 


and 


(5) 


r FF x 


tsy \ 

(F m -£,)]+ xli + SIF Q {h '%,+ SIF Q ,F ’M i - sll Q iF \lk F ) =Q, (Eq. 193) 

V i 


where [cf. Eq. (171)] 


(5) ,V 


(Eq. 194) 


These equations must now be placed into a usable state-space form. 


Define the following: 


=1 


<7i 

Qi 

<73 


(Eq. 195) 
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(Eq. 196a, b, c) 


and 



= -(O h x c + (0 h {i<l F ~L FEy -l i Q^\ (Eq. 197a, b, c, d) 



(Eq. 198) 


(Eq. 199) 


u = I. 


(Eq. 200) 


Note that Equations (196) and (199) express Kane’s Kinematical Equations. 

Using state-definition equations (195), (196), and (199); force equations (177) and (180); 
acceleration equation (192); and disturbance equations (40) and (41); Kane’s Dynamical Equation 
(191) becomes 


P^kt 4, »' 5 ^»yV+ h, 

m f m m l m 


/it \ \ 

(S U 4 J S) n -I- 1 (S) F 4 - L F (S) r 
- Q.in+ Qld+ — Lb+ — b uta Lb- 

m m 


Likewise, using state-definition equations (195), (196), (198), and (199); force equations (177) and 
(180); moment equations (82), (83), (181), and (184); angular acceleration equation (194); and 
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disturbance equation (42); Kane’s Dynamical Equation (193) becomes 

i,= ( s ' r Qr‘M„)x . +( m Qr l Mjx„ +( ~x_ a 

+ ( S ' F G/ _1 M„,- S ' F G/ _ 'M,/ /F e '"C-f. S ' r Q T )i,+ ( Slr QI~'M c )u Eq. (202) 

+ ( S "Q /-M J“v. + ( s,f Q )( f, «, + /" '"«,)■ 

In terms of Euler parameters x d ~ {s r 'P [where 7 s p is defined by Eq. (26)], for small rotation 

angle <|) and negligible angular velocity N co s the rotational states are 

x_ d = 2x d , (Eq. 203) 

and x_ e = 2x e . (Eq. 204) 

Kinematical equations (196) and (199), and dynamical equations (201) and (202), now reduce to 
the forms found previously, viz.. Equations (123), (126), (124) and (127), respectively. 

Concluding Remarks 

This paper has presented the derivation of algebraic, state-space equations for the Canadian 
Space Agency’s Microgravity Vibration Isolation Mount. The states employed include payload 
relative translational position (x ^ ) and velocity (xf), payload relative rotation (x d ) and rotation 

rate (jff), and payload translational acceleration (vf ) . Feedback of x^ : corresponds to a change in 
effective umbilical translational stiffness, with the effective umbilical assumed to be attached at the 
flotor center of mass. Similarly, feedback of x£ , x d , or x_ c e corresponds, respectively, to a change 

in translational damping, rotational stiffness, or rotational damping, for the same effective 
umbilical. Likewise, feedback of payload translational acceleration causes a change in effective 
payload mass. Thus, a cost functional which penalizes these states produces an intuitive effect on 
system effective stiffness, damping, and inertia values. 
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The acceleration states can be selected to pertain to any arbitrary point on the flotor. This 
allows an optimal controller to be developed which penalizes directly the acceleration of any 
significant point of interest, such as the location of a crystal in a crystal-growth experiment. 

The equations have been put into state-space form so that the powerful controller-design 
methods of optimal control theory (e.g., H 2 synthesis, H x synthesis, jJ synthesis, mixed-/! 
synthesis, and p analysis) can be used. References [1 1], [12], and [13] detail the H 2 optimal 
controller design approach used for MIM. 
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